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h^ ; Abstract 

\Q ' We evaluate the rare radiative kaon decays 1£l,s ~~ ► ^7- Applying the requirements of gauge invariance 

and Bose symmetry, we derive a general form of the decay amplitude, including both parity-conserving 

i—i. and parity- violating contributions. We employ a chiral-Lagrangian approach combined with dimensional 
analysis arguments to estimate the branching ratios of these decays in the standard model, obtaining 
values as large as B(K L — > 3j) ~lx 10 -14 and B(K S — > 37) ~ 2 x 10~ 17 , which exceed those found 
previously by a few orders of magnitude. Measurements on the branching ratios which are significantly 

r^, ■ larger than these numbers would likely hint at the presence of new physics beyond the standard model. 
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I. INTRODUCTION 

The rare kaon decays into three photons, K^ — > 37 and K$ — > 37, can occur in the absence of 
CP violation. One might then naively expect from the measured branching ratio B{K\^ — > 27) ~ 
5.5 x 10" 4 [1] that B(K L -> 37) ~ a cm B(K L -> 27) ~ 4 x 10 -6 . However, this expectation is 
already too large in comparison to the result of the first experimental search for the 37 mode, 
B{K\, — > 37) < 2.4 x 10~ 7 [1, 2]. As for Kg — ¥ 37, there is currently no experimental information 
available on it, but it is likely to be more suppressed than expected as well. 

It turns out that the considerable smallness of the K — > 37 rates has to do with the constraints 
imposed on the decay amplitude by gauge invariance and Bose symmetry [3]. Gauge invariance 
implies that the total angular momentum J of any two of the three photons in the 37 final-state 
cannot be zero, whereas Bose statistics forbids the photon pair to have J = 1. Since each of 
the photon pairs must have J > 2, the decay amplitude suffers from a large number of angular 
momentum suppression factors. 

The rates of i^L,s —* %1 were roughly estimated a while ago in Ref. [3]. The calculation 
was based on a simple model in which K — > 37 is assumed to proceed from K — > 7r°7r°7 with 
the 7T° pair immediately converting into a photon pair. The resulting branching ratios are tiny, 
B(K L -> 3 7 ) ~ 3 x 1(T 19 and B(K S -»■ 3 7 ) ~ 5 x 1(T 22 [3]. 

Here we take another look at these rare decays, partly motivated by a new search for K^ — > 37 
currently being performed in the E391a experiment at KEK [4]. Since the estimate obtained in 
Ref. [3] resulted from only one contributing diagram, it is possible that other contributions exist 
which can enhance the decay rate. In this study we start with a general form of the K — > 37 
amplitude subject to the restrictions from gauge invariance and Bose statistics. We then use 
a chiral-Lagrangian framework along with dimensional-analysis arguments to explore the size of 
various important contributions. This finally leads us to i^L,s — > 37 rates which exceed those 
previously estimated by a few orders of magnitude. 

II. DECAY AMPLITUDES AND RATES 

The decay K — > 37 being a weak transition, its amplitude hA[K — > 37) generally consists of 
separate terms describing the parity conserving (PC) and parity violating (PV) components of 
the transition. Accordingly, one can write Ai (K — > 37) as a sum of their respective contributions, 

M{K -)■ 3 7 ) = M£ c + Mgy , (1) 

KA K — p* p* p* A/f a/3M KA K — p* p* p* M a/5fl (9) 

where e l 2 3 are the polarization vectors of the photons. 

To construct a general form of the K — > 37 amplitude, one follows a well-known prescription. 
The presence of photons implies that M.(K — ¥ 37) has to be gauge invariant. With three photons 
in the final state, Bose statistics dictates that the amplitude be symmetric under interchange of 
any two of the photons. These requirements must be satisfied by .Mpc and A^pv separately. 



We deal with Mpy first, as it turns out to be simpler than Mpc and is also more relevant 
to Kl — y 37, which is our main decay of interest, relegating some details to Appendix A. We 
also assume that the photons are on-shell. Thus we get 

Mpf = (g^z - k%k?){k»y - k^ x) F(x,y,z) 
+ {gtoy - kik£) (A£ x - k% z) F(z, x, y) 
+ (g a »x - fc 3 X) (k§ z-k{y) F(y, z, x) 
+ [g^ (K y - e 2 x) + gf* (k« x - fc« z) 

+ g a » (kl z-kfy) + k%k(k% - *£*£*£] G(x, y, z) , (3) 

where k 1 2 3 are the momenta of the photons with polarizations e 1 2 3 , respectively, 

x = k 1 ■ k 3 , y = k 2 ■ k 3 , z = k 1 ■ k 2 , (4) 

and the functions F and G must be free of kinematic singularities and satisfy the relations 

F(u,v,w) = —F(v,u,w), 
G(u,v,w) = —G(v,u,w) = —G(w,v,u) = —G(u,w,v) , (5) 

with u,v,w each being any one of the invariants k i ■ k». This amplitude agrees with the one 
derived in Ref. [5] for 7r° — y 37. 

For the parity-conserving contribution, the form with the desired symmetry properties can be 
expressed as 

Mpf = [{g aP z - k«k{) e^ T F(x, y, z) + (g^y - fcf A£) e a ^ F(z, x, y) 
+ (g^x - kSk?) e^ T F(y, z, x)] k lp k 2a k 3r 
+ [(e 2 k{ - k?k T 2 ) t a ^ H(x, y, z) + (kfJ4 ~ k%k p 3 ) e^ T H(z, x, y) 
+ [k l k 3 — k 3 k 1 ) e rL{y, z, x)\ k Xp k 2a k 3r 

_|_ 1 I ap e iip<7T _|_ pcr^a/ipr _|_ fip^apor _ aa^ppr 

1 qPv^apcTT _|_ ar^otPpp 1 pa^aBpr _ Br^appa 
+ g ap e /3par + g pr ^Bpa + g pp ^ar _ g ar ^pvj k lp k 2ff k 3r G(x, y, Z) , (6) 

where the functions F, Q, and % are also free of kinematic singularities and satisfy 

F(u,v,w) = —F(v,u,w) , %(u,v,w) = —%(v,u,w) , 

G{u,v,w) = —Q(v,u,w) = —Q(w,v,u) = —Q(u,w,v) . (7) 

The formula for Mpc above may have been constructed for the first time in this paper. 

After summing |A4pV + A4p C | 2 over the photon polarizations, we find that there is no inter- 
ference between the PC and PV contributions in the result, in accord with expectation. It is 
given by 



\M(K^3r/)\ 2 = ^(1-Mfvl 2 + \M&\* 



E 

pol pol 



where 

EI-^pvI 2 = 4{l*ifV + I^IV + l^l 2 ^ 2 + 2 |Gf (9) 

po1 + Re [F*F 2 yz + F;F 3 x y + F*F ± xz + 2 (F*z + F*y + F*x) G] }x y z , 

El-^pcl 2 = 4 ( (l^il 2 + l«il V + (l^l 2 + I« 2 IV + (l^l 2 + |« s l V + 2 \gf (io) 

po1 + Re[(^ + HI) {T 2 + H 2 + 2g/y)yz+ (J> + H* 2 ) (T 3 + H 3 + 2Q/x)xy 

F x = F(x,y,z) , F 2 = F(z,x,y) , F 3 = F(y,z,z) , (11) 

similarly for ^23 and H x 2 3 , and 

G = G(x,y,z), G = G(x,y,z) . (12) 

The resulting decay rate can be written as 

r(X->3 7 ) = I j, f ds 12 ds 2 ^\M{K ^ ^W (13) 

K ' J pol 

where the 3! accounts for the three photons being identical particles and s mn = (k m + k n ) 2 . Using 
the formulas above, we provide our numerical estimates in the next section. 

Before moving on, we remark that the expressions in Eqs. (3), (6), (9), and (10) apply more 
generally to the decay of other neutral pseudoscalar particle into three photons. Furthermore, 
they also work for a neutral scalar particle decaying into three photons, but with the PC and PV 
contributions interchanged. 

III. ESTIMATE OF K L , S -*■ $1 RATES 

To explore the size of the leading contributions to M.(K — > 37), we adopt a chiral-Lagrangian 
approach [6]. Accordingly, they are expected to arise from the relevant terms in the chiral ex- 
pansion and give rise to terms in the functions F, G, J 7 , g, and H with the lowest numbers of 
powers of the photon momenta k i . Since there are in principle many possible contributions to 
the amplitude, from tree and loop diagrams, with mostly unknown parameters, we consider a few 
representative contributions and rely on dimensional- analysis arguments to determine their size. 

A. Kl -»■ 37 

Neglecting CP violation, we can concentrate on the PV part of the amplitude, in Eq. (3), 
as i^L — > 37 violates charge-conjugation invariance. Thus, for F and G satisfying Eq. (5) we find 
the simplest form 

F(u,v,w) = c F (u — v) , 
G(u, v, w) = c G [(u - v)f(w) + (v- w)f(u) + {w- u)f(v)] , (14) 



where c FG are constants and / is any well-behaved function, although it cannot be a constant 
if G is to be nonzero. This implies that F and G contain at least two and four powers of k t , 
respectively, as the momentum power in the chiral expansion in the meson sector always increases 
by even numbers. It follows that Mpy in Eq. (3) contains at least seven powers of k { . 

To assess the leading contributions to Mpy, we first consider an example of a weak chiral 
Lagrangian for strangeness changing, \AS\ = 1, transitions within the standard model (SM) 
which is odd under parity, has seven derivatives, and couples a kaon to three photons in a gauge- 
invariant way. As is well known, the weak chiral Lagrangian for such transitions in the SM is 
dominated by contributions which transform as (8l, 1r) [6] and has to be invariant under the CPS 
transformation [7], which is the product of the ordinary CP transformation and the switching of 
the s and d quarks. An example with the required properties is 

£pv = o, ft^ (VMT) [W' V^ + (V^JW] W p ) + H.c. 

= ^YTf- 9 a F^ {d a F pa + d a F pa ) d»K° d°F pv + ■■■ + H.c. , (15) 

where c 7 is a constant, F pu = d /Ji A l/ — d v A p is the usual photon field strength tensor, and only the 
relevant part is displayed in the second line. In the first line, £ = e'*''' 2 ^ contains the lightest octet 
of pseudoscalar mesons via the 3x3 matrix tp [6], with f = f n = 92MeV being the pion decay 
constant, h is a 3x3 matrix having elements h kl = S k2 S 3l which selects out s — > d transitions, 
Vnv = e {t^Qi + £Q^)Fij,v with Q = diag(2, — 1, — 1)/3 being the quark charge matrix, V ' a X = 
d a X+±[&d a Z+Zd a &+ie(ZQ?+?QZ)4*,X], fmdU a = i(?d a S-Zd a ?)+e(ZQ&-?QZ)A a . 
The Lagrangian in Eq. (15) yields 



32y / 2 ic e 3 
F(u,v,w) = — — j 1 — (u-v), G(u,v,w) = (16) 

in A4(K°) and the same functions in .A/f (-ft' ), but with c 7 replaced with c 7 if CP violation is 
present. 

If CP is conserved, only .Mpy contributes to the K L — >■ 37 amplitude. In that case, upon 
applying Eq. (16) in Eq. (9) and adopting the convention K L = (K° + K } /\/2, we find 

El A A I TS > o ^l 2 C 7 e /22, 22, 22 2 2 2 \ ;,m 

\M{K L -^37)1 = , 2 [x y +y z +x z -xyz -xy z-x yzjxyz , (17) 

poi il ^^ 

having made use of Eq. (4). To calculate the decay rate, we then need the value of c 7 . Since it is 
not possible at present to determine this constant rigorously from the quark-level parameters, we 
estimate it with the aid of naive dimensional analysis [8]. Thus we obtain the order-of- magnitude 
value 

^ G F \ C ti „ 1Q 10 -9 Ge y-6 (lg) 

where A c = 0.22 is the Cabibbo mixing parameter and A represents the scale at which the 
chiral Lagrangian approach breaks down, and so we take A = m = 775 MeV [1]. The resulting 
branching ratio is 

B(K L ->■ 3 7 ) ~ 7.4 x 10" 17 . (19) 



This number is about 250 times larger than the earlier prediction, B(K^ — > 37) ~ 3 x 10~ 19 , 
given in Ref. [3]. We have repeated that calculation and come to a different result, which we 
briefly discuss here. In the model employed in Ref. [3], the amplitude Ai(K L — > 37) is roughly 
approximated by 



/l 



X 



L 



m^ 



(£3 'Pi ^3"PX - ^'Px *VPl) *V ( 2 Pl - P«") 



G s 12 1 fej/Pi k 2 ■p 1 ki'Pi 1 ^2' Pi 1 I *// */' 

~^2~ F~t ^ + PlMPli' TT1T K 2^Plu -j—jT K luPln I £ i e i 

" L p \ h l ^2 h l ^2 "a ^2 

+ [permutations of (e 1; A^), (e 2 ,k 2 ), (£3,^3)] , (20) 

coming from a 7T°-loop diagram, where p K is the K^ momentum, p 1 2 are the 7T° momenta in the 
loop, h L ~ 6.6 x 10~ 8 , and G = y g£ with g W7r7 = 0.77 GeV -1 . To compute this requires the 
evaluation of the integrals 

K^ = [ P±P± *< 4 > (P- Pl - p 2 ) ptpWPlf{Pi ■ ft) , 

J A PlO z fto 

L^ = [ pl d ^l *W {P- Pl - p 2 ) Atfj&Zffa ■ p 2 ) . (21) 

J A PlO Z fto 

We have collected the results in Appendix B, where the expression for K^ vpu agrees with that given 
in Ref. [3] , but our L^ up<T differs from theirs. With Im Ai (K L — > 37) computed and subsequently 
equated to Aipy(K L ) = £ ia e 2/3 e 3uMpy , we then extract 



F(u,v,w) 
G(u,v,w) 



^-( v -2mlf /2 e(v-2ml) - -±= (u - 2m\f 2 6(u - 2m;) 
'v \/U 



"' ~ V (u-2mlf 2 9(u-2ml) + ^^( v - 2mlf 2 9(v - 2ml) 



u 



V — U, _ 2 x5/2 



+ "—l( w -2miy z e(w-2ml) 

\/W 



(22) 



Glh 6.3 x 10- 9 GeV- 9 . (23) 



120 71771^171^ 



Incorporating these into Eqs. (9) and (13) leads us to B[K\, — ¥ 37) ~ 1.0 x 10 18 , which is three 
times greater than the number found in Ref. [3]. 

Since the results in the preceding paragraph arise from a loop diagram [3], it is instructive to 
compare them to tree-level contributions of the same chiral order. From the F and G formulas in 
Eq. (22) we can see that they contain four and six powers of the photon momenta k i , respectively, 
and hence their contributions to M (K L — ¥ 37) have nine powers of k i . If a weak chiral Lagrangian 
with nine derivatives contributes to this amplitude at tree level, the size of its coupling constant 
in F(u,v,w) is predicted by naive dimensional analysis to be 

_ e*G ¥ \ c fl ^ 10 9 

v^A 10 



This is about 12 times smaller than c in Eq. (23). One may then suggest based on this comparison 
that loop contributions with seven powers of k i could also be enhanced by a similar amount relative 
to the tree-level contribution from Lagrangians with seven derivatives, such as that in Eq. (15). 
If such enhancement occurs, we may have 

B(K L -> 3 7 ) ~ 1 x 1(T 14 , (25) 

instead of Eq. (19). The numbers in Eqs. (19) and (25) can be taken to be representative values 
coming from the various contributions to this decay and also to indicate the level of uncertainty 
involved in our crude calculation. It is therefore reasonable to conclude that 

7 x 1(T 17 < B(K L -> 3 7 ) < 1 x 1(T 14 . (26) 

If CP violation is not neglected, then .Mff c also contributes to the K h — > 3 7 rate, via the 
e K term in K L ~ [K° + K° + e K (K° - K°)]/V2, where e K ~2x 10" 3 is the CP- violation 
parameter in kaon mixing. Now, the similarity between Eqs. (5) and (7) suggests that F and % 
(Q) are comparable in size to F (G). We can then expect that -M^ c is, at most, also comparable 
to -Mfy- This implies that the effect of CP violation on B[K L — > 3 7 ) is small, and hence the 
predicted branching ratio is still what is quoted in Eq. (26). 



B. K s -*■ 3 7 

The rate of this decay is determined mostly by the parity-conserving contribution A4p C . The 
branching ratio predicted in Ref. [3] is B(K$ — > 3 7 ) ~ 5 x 10~ 22 . Repeating the calculation, as in 
the Kl — > 3 7 case, we find instead a value three times larger, B^Ka, — > 3 7 ) ~ 1.8 x 10~ 21 . Since 
this arises from a loop contribution involving nine powers of the photon momenta k^ we need to 
consider as before the lower-order contributions, with seven powers of k t) which may be larger. 

We take the leading-order form F{u,v,w) ~ 7i(u,v,w) = c(u — v), satisfying Eq. (7), with c 
being a constant and Q = 0. In this case, the situation is similar to that of M.?c w hh F and G 
given in Eq. (16). More precisely, making a comparison of Z) po il-^Pcl 2 m Eq- (10) and I] po i|A^p V | 2 
in Eq. (9) for the two cases, respectively, one can see that their decay distributions have the same 
functional dependence on x, y, and z. It follows that T(K S — > 3 7 ) can be expected to be roughly 
of the same order as T(K L — > 3 7 ). Interestingly, the measured rates of the corresponding 2 7 
modes are also of similar order, r(-K" s — > 2 7 ) ~ 2.7 T(K L — > 2 7 ) [1]. In view of Eq. (26), we can 
therefore predict that 

1 x 1(T 19 < B(K S ->■ 3 7 ) < 2 x 10" 17 . (27) 

IV. CONCLUSIONS 

We have revisited the rare kaon decay K — > 3 7 , which is expected to be much suppressed 
because its amplitude has a large number of angular momentum suppression factors. We have 



constructed a general form of the decay amplitude which satisfies the requirements of gauge invari- 
ance and Bose symmetry and includes both parity-conserving and parity-violating contributions. 
We have in addition calculated the squared amplitude, summed over the photon polarizations, 
which can be useful to produce a Dalitz plot distribution of the decay. These results are applica- 
ble generally to the decay of any spinless particle into three photons. 

More specifically, we have dealt mainly with K^ — y 37, which is currently the subject of a new 
experimental search at KEK, but also evaluated Kg — * 37, albeit more briefly. To explore the 
leading-order contributions to their amplitudes, we have adopted a chiral-Lagrangian approach in 
the context of the standard model. This implies that there are many possible contributions to the 
amplitudes, from tree and loop diagrams, with mostly unknown parameters. Consequently, for 
definiteness we have considered a number of representative contributions and used dimensional- 
analysis arguments to estimate their size. This has finally led us to arrive at branching ratios as 
large as B(K L — y 37) ~lx 10~ 14 and B(K S — > 37) ~2x 10~ 17 , which exceed those estimated 
before by a few orders of magnitude, but are still very small. Nevertheless, any experimental 
findings on the branching ratios which are significantly greater than these numbers would likely 
signal the effect of new physics beyond the standard model. 
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Appendix A: Derivation of amplitudes 

To obtain a general form of the K — y 37 amplitude with the desired properties, we start with 
expressions for Mpy and Mp^ consisting of terms involving all possible combinations of the 
available tensors, A;" 2 3 , g VK , and ^ V(JT . Thus we have 

Mpf = g afi {K a, + fc£ o 2 ) + g afl (fcf a 3 + fc£ a 4 ) + <A (k° a 5 + k$ a 6 ) (Al) 

+ k«[k{ {k{\ + k^b 2 ) + k§(k?b 3 + k%b A )} + k«[kP(k»b 5 + k^b 6 ) + 4{ktb 7 + k»b 8 )} , 

MgT = e a ^ (k lpCl + k 2p c 2 + k 3p c 3 ) + (g^e^d, + g^e^d 2 + g^e a ^d 3 )k lp k 2a k 3r 

+ [{k a 2 A + k« f 2 y«° + (k( f 3 + kg f A )e*» + (a? u + K hy aPp<T ] K k 2« 

+ [ (A£ 9! + K g 2 ) #"» + {K g 3 + k» g 4 ) e^ + (k?g 5 + k§ g 6 ) <™"] k lp k 3a 

+ [(^ h x + k» h 2 )e a ^ + (k§ h 3 + k[ h 4 y^ + (ks h + k« \y^\ k 2p k 3a 

+ J yk 2 yk 1 1 1 + k 3 1 2 j + k 3 yk 1 1 3 + k 3 / 4 J J e + yk 2 [k 1 l 5 + k 2 / 6 J + k 3 [k 1 L 7 + k 2 1 8 ) J e 
+ [/c 1 [^Ig + k 2 l w ) + k 3 (k 1 l 11 + k 2 l 12 )\e \k lp k 2a k 3T , (A2) 

where the a { , b iy q, d t1 f iy g iy h t , and l i are functions dependent on the invariants k i • k-. After 
the requirements of gauge invariance and Bose symmetry have been imposed on Mpy and Mp^ 



separately, they each contain a much smaller number of functions. For M^q , we make further 
simplification with the aid of Schouten's identity, which states that in four dimensions a tensor 
with five or more Lorentz indices vanishes identically if it is completely antisymmetric with respect 
to five or more of the indices. 1 Such a tensor is 

ap^upar _ au^ppar _ ap^upar _ aa^i/pfxr _ ar^upap, _ q (A3) 

which is fully antisymmetric with respect to fx, u, p, a, r. We display our results for Mpy and Mpc 
in Eqs. (3) and (6), respectively, in the case of on-shell photons. The Q terms in Eq. (6) could be 
simplified further using Eq. (A3), but then they would not be manifestly Bose-symmetric. 

Appendix B: Integrals 

The integrals in Eq. (21) can be written in terms of all the possible appropriate combinations 
of the available tensors, g al3 and P v , as 

' A 5 (s,mf,m|) 



71 . 
J^pvpa _ 



10 S 



{g pu g pa + g^g va + g fia g' /p )K 1 + ppp»ppp° k 3 

+ {g^PPP* + gWp v P° + g^P^P + gVPp^pv ( B1 ) 

+ g^P^pp + g^ppp") K 2 ] f(l{ s -ml- m 2 .)) , 



L ^ P c = K}H s > m i m i) 



10 s 



(gWgP* + gWg™ + g v Pg^)h x + pP>P»pPp° L 4 

+ (g^PPP* + gwp v p° + g v Pp v P°) L 2 (B2) 

+ {f°P v pi> + g^pppp + gP°pPp v )Lh(\(s - m\ - m 2 .)) , 

where X(u,v,w) = u 2 + v 2 + w 2 — 2uv — 2v w — 1u w, m 2 = p 2 , m 2 = p 2 , and the coefficients 
K 1 2 3 and L 1 2 3 4 are functions of s = P 2 and m 1 2 . We then derive 

' [s 2 -2s(m 2 + m 2 ) + (m 2 -m 2 ) 2 ] 2 , 

[s — {m 1 + m 2 ) 2 ] [s — {m 1 — m 2 ) 2 ] [3s 2 + s (4m 2 — 6m 2 .) + 3(m 2 — "ml) 2 ] , 

{s 4 + s 2 (m\ — Qm\m\ + 6m 2 ) + s(m 2 — 4m 2 ,) [s 2 + (m 2 — "ml) 2 ] + {jn\ — m 2;) 4 } 5 

(B3) 

[s 2 - 2s (m 2 + m 2 ) + (m 2 - m 2 ) 2 ] 2 , 

[s — (m 1 + m 2 ) 2 ] [s — (m x — m 2 ) 2 ] [2s 2 + s (m 2 + m 2 ) — 3(m 2 — m|) 2 ] , 

1 Some other examples of the use of Schouten's identity can be found in Ref. [9] . 



jx 1 - 


48 s 2 


K 2 = 


= - 1 

24 s 3 


K 3 = 


s 4 l 


Li - 


-1 

48 s 2 


L 2 - 


-1 
24 s 3 



L 3 = 7^-^ [s - (toj + m 2 ) 2 ] [s - [m l - m 2 ) 2 ] [3s 2 + s (Ami ~ 6 m i) + 3(m 2 - m 2 ,) 2 ] , 

L 4 = — - [s 4 + s 3 (m 2 + mf) + s 2 (m\ + Am\m\ - 9m2) 

+ s (ml + 1 lm 2 ) (m\ - m 2 ) 2 - 4 (m 2 - m 2 ) 4 ] . (B4) 

Our formula above for K^ vpa (L^ upa ) agrees (disagrees) with that found in Ref. [3]. In obtaining 

Eq. (22), we set m 1 = m 2 = m n . 
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